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Feshbach resonances in ultracold ^^Rb— ®^Rb and ^Li— ^^Rb mixtures 
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We present an analysis of experimentally accessible magnetic Feshbach resonances in ultracold 
heteronuclear *^Rb-*^Rb and ^Li-*^Rb mixtures. Using recent experimental measurements of the 
triplet scattering lengths for ^Li-*^Rb and ^Li-*^Rb mixtures and Feshbach resonances for one 
combination of atomic states, we create model potential curves and fine tune them to reproduce the 
measured resonances and to predict the location of several experimentally relevant resonances in 
Li-Rb collisions. To model *^Rb-*^Rb collisions, we use accurate Rb2 potentials obtained previously 
from the analysis of experiments on *^Rb-*^Rb collisions. We find resonances that occur at very 
low magnetic fields, below 10 G, which may be useful for entanglement generation in optical lattices 
or atom chip magnetic traps. 

PACS numbers: 34.50.Cx, 34.20.-b 
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I. INTRODUCTION 

The field of ultracold atomic gases has been revolu- 
tionized by the discovery and use of magnetically tun- 
able Feshbach resonances jT], [3, These resonances 
occur when the energy of a quasibound molecular state 
( "closed" channel) is tuned into degeneracy with the en- 
ergy of a colliding atomic pair in an "open" channel by 
an externally applied magnetic field. Near a resonance, 
the coupling between the states strongly affects both the 
elastic and inelastic collision cross sections of the colliding 
atoms. First observed at low temperatures in homonu- 
clear collisions of Na atoms [41] , Feshbach resonances pro- 
vide a way to tune the microscopic interactions and to co- 
herently (reversibly) form ultracold molecules These 
features have been exploited for the study of a wide vari- 
ety of phenomena in ultracold atomic gases including the 
controlled collapse of a Bose-Einstein condensate (BEC) 
0, bright solitons 0, Q, and BCS superfluid pairing of 
a Fermi gas 

Mixtures of different ultracold atoms are interesting 
objects of study as they may exhibit scattering proper- 
ties different from those of indistinguishable species. For 
example, broken exchange symmetry allows for s-wave 
scattering of different fermions in the same hyperfine 
state. In addition, the Feshbach spectrum of heteronu- 
clear molecules is more complex because the hyperfine 
constants are unequal for the two atoms and there are 
four rather than three zero-field collision asymptotes [l^. 
Heteronuclear Feshbach resonances have recently been 
observed in *^U-^^^a [H, ^^K-^'^Rh [3, ^SRb-S^Rb 
[H, 6Li-87Rb PI, and 6Li-40K collisions These 
resonances provide a mechanism to tune interspecies in- 
teractions and may enable the study of boson-mediated 
cooper pairing 18] and supersolid order |19J predicted 
to occur in Bose-Fermi mixtures. Feshbach resonances 
in heteronuclear systems can be used to produce ultra- 
cold polar molecules [l^] . The dynamics of ultracold po- 
lar molecules is determined by tunable and anisotropic 
electric dipole-dipole interactions, and the creation of ul- 



tracold polar molecules is predicted to allow for major 
advance in several different fields of physics and chem- 
istry. 

In this paper, we analyze collisions in ultracold *^Rb- 
^^Rb and ^Li-^^Rb mixtures. Recently, Tiesinga com- 
bined the results of spectroscopic and scattering exper- 
iments on ®^Rb ensembles to produce a new set of ac- 
curate Rb2 interaction potentials of singlet and triplet 
symmetry [21]. With the new potentials, Tiesinga was 
able to reproduce the experimentally observed Feshbach 
resonances in ultracold collisions of ^"^Rb atoms with a 
high level of accuracy. Here, we extend this work to 
heteronuclear rubidium dimers. We calculate the mag- 
netic field dependence of *^Rb-*^Rb elastic scattering 
cross sections using the new interaction potentials and 
find good agreement with the recent experiment of Papp 
and Wieman [l^ and the theoretical analysis of Burke 
et al. 0. We also study ultracold scattering of Rb 
isotopes in hyperfine states other than the fully spin- 
stretched state |2,— 2)85 ® If;— 1)87 and find resonances 
that occur at low magnetic fields, below 10 G. For the 
Li-Rb system, we create model potentials combining the 
results from recent ah initio calculations for the LiRb 
dimer [2^, and then fine tune them to reproduce the 
experimentally measured resonances. For this purpose 
we utilize an asymptotic bound state model [l3, 17, 2j| 
to determine the energies of the last bound states of the 
triplet and singlet potentials consistent with the experi- 
mentally observed location of two Feshbach resonances. 
Guided by the results of the model analysis, we gener- 
ate the corresponding potentials and compute both the 
location and widths of the Feshbach resonances using a 
full quantum scattering calculation. This process allows 
us to systematically converge to the optimal singlet and 
triplet potentials for the Li-Rb dimer. These potentials 
fully characterize the scattering properties in any combi- 
nation of atomic spin states and allow us to predict both 
the triplet scattering length and the location of several 
experimentally accessible resonances for this mixture. 
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A. Calculation method 

The Hamiltonian for the collision of two alkali metal 
atoms is H = Hj-^i + Vm + Vd + Vb where the first term 
accounts for the relative motion of the atoms, Vjif models 
the hyperfine interactions, and Vb models the interaction 
of the collision complex with the external magnetic field. 
We have verified that the magnetic dipole interaction, Vd, 
has no effect on the observables described in this paper 
so we neglect it. 
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where fi is the reduced mass of the atoms, R is the inter- 
atomic distance and / is the rotational angular momen- 
tum of the collision complex. The total scattering wave 
function is expanded in a fully uncoupled, space-fixed 
basis set: 
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where Fa.i.Mi{R)\lMi) are the radial basis states and 
|a) \IaMO\SaMs^)\IbMi^)\SbMs^) are the atomic 
spin states. The electron and nuclear spins for atom 
a (b) are denoted by Sq (Sb) and la (Ifc). The to- 
tal electron spin for the complex is the vector sum of 
the electron spins of atoms a and 6, S = + S^. 
V{R) — J2s^s(^)-^s describes the inter-atomic inter- 
action potential where Ps = J2ms \SMs){SMs\ is the 
projection operator onto the singlet S' = or triplet 
S = 1 electron spin configuration of the molecule and 
Vs (R) is the corresponding interatomic potential for the 
singlet X^T, or triplet a^S state of LiRb. The interaction 
of the collision complex with the external magnetic field 
is given by 

Vb - 5e(Sa + Sb) . B ~ {gi^K + .gi'^Ifa) ■ B (3) 

where ge and gn are the electron and nuclear gyromag- 
netic ratios. The hyperfine interaction is described by 



(4) 



where 7a and jb are the hyperfine interaction constants 
for atoms a and b. 

The matrix elements of the Hamiltonian in the fully 
uncoupled basis are determined as described in Ref . [2^ , 
and the resulting coupled-channel equations are prop- 
agated into the asymptotic region where the standard 
asymptotic boundary conditions [l^l are applied to ex- 
tract the scattering matrix elements. The integral cross 
sections and the scattering lengths are computed from 
the ^-matrix elements as described in Ref. [28| . 

Using this procedure, we calculate the elastic scatter- 
ing cross sections as functions of the magnetic field at a 
fixed collision energy. To ensure the convergence of the 
calculations, we employ a dense propagation grid from 2.0 



to 800 OB (flB = 0.0529177 nm) with a step of 0.005 or. 
The s-wave scattering length near a Feshbach resonance 
has the form (H 



a{B) = flbg 1 
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where Bq is the position and AB is the width of the 
resonance. The background scattering length Obg does 
not depend on the magnetic field. For each resonance, 
we extract the parameters Bq and AB from the mag- 
netic field dependence of the scattering length using Eq. 
[5] In the case of the ^^Rb-^^Rb resonances, we have ver- 
ified our results for Bq by analyzing the derivative of the 
eigenphase sum [2^. The resonance positions obtained 
in this way are found to be in excellent agreement with 
the results of Eq. [5l 



II. LOW FIELD RESONANCES IN **^RB- 
MIXTURES 
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In 2006, Papp and Wieman created heteronuclear 
molecules in a mixture of ^^Rb and ^^Rb atoms by a 
linear ramp of the magnetic field and by resonant-field 
modulation [15j. For this purpose they used two Fes- 
hbach resonances in the |2, — 2)35 ^ |1,— 1)87 magnetic 
hyperfine state (where |/, my) is the usual notation for 
the atomic hyperfine states). Burke et al. considered 
magnetic Feshbach resonances in the weak magnetic field 
seeking states of *^Rb-^^Rb and suggested optimal con- 
ditions for sympathetic cooling of ®^Rb atoms by elas- 
tic collisions with evaporatively cooled *^Rb [l^l- Using 
accurate Rb2 potentials of singlet and triplet symmetry 
[2l|, we analyze the magnetic field dependence of *^Rb- 
^^Rb elastic scattering cross sections. Our primary moti- 
vation for this work was to search for low field Feshbach 
resonances which would be candidates for the realization 
of RF induced Feshbach resonances using large RF or 
microwave fields generated near the surface of an atom 
chip ^30| . 



A. Results 

Figure [T] shows the s-wave elastic scattering cross sec- 
tions for *^Rb-*^Rb collisions in the weak magnetic field 
seeking state |2, — 2)85 ® |1,— 1)87- Two resonances are 
immediately apparent from this graph: the one at lower 
energy (resonance I) and a broader peak around 380 G 
(resonance II). These resonances are of significant inter- 
est for sympathetic cooling of ^^Rb by collisions with 
ultracold ^^Rb [l^, [l^], and it is important to know 
their paramaters with high accuracy. Table |T] compares 
our calculated resonance positions and widths with re- 
cent experimental results of Papp and Wieman [Tsj . The 
agreement is generally good, although the calculated po- 
sition of resonance I (271.05 G) is shifted to higher mag- 
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Magnetic field (G) Magnetic field (G) 

Figure 1: Elastic scattering cross sections for s-wave colli- 
sions of *^Rb and *^Rb atoms in the lowest weak magnetic 
field seeking hyperfine state |2, — 2)g5 ® |1,— 1)37 versus the 
applied magnetic field (upper panel). The collision energy is 
144 nK (see text). The lower panels show more details of the 
resonances. 
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Figure 2: Magnetic field dependence of s-wave elastic scat- 
tering cross sections for *®Rb-*'^Rb in different hyperfine 
states. Full line - |2,0)g5 ® |1,— l)g7, long dashed line - 
|2, -1)85 ® |1, -1)87, short dashed line - [2, 1>85 ® |1, -1)87- 
The collision energy is 144 uK. 



TABLE I: Positions and widths of magnetic Feshbach res- 
onances in ^^Rb-^'^Rb collisions for different initial atomic 
hyperfine states. All atomic states considered are from the 
lowest hyperfine manifold (/ = 1 for ^'^Rb, / = 2 for *^Rb). 
The numbers in parentheses are taken from calculations of 
Burke et al. [H. 
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netic fields by 5.6 G with respect to the experimen- 
tal value The calculated width of resonance I is 
AB = 3.66 G, consistent with the experimentally mea- 
sured width, AScxp = 5.8 ± 0.4 G, of the trap loss fea- 
ture [15]. Although they are related, the width of the 
trap loss feature, ABcxp, is not equivalent to the width 
of the scattering length singularity AB Table U also 
shows the results obtained with a different set of inter- 
action potentials by Burke et al. [22]. The lower bound 
of their estimate almost coincides with the experimental 
prediction, whereas the upper bound is only 2 G away 
from our result of 271.05 G. 

The second resonance in Fig. [1] is a broad resonance 
centered at 382.7 G with a calculated width of 33.2 G. 
Table U shows that our calculations overestimate the res- 
onance position by 10 G. On the other hand, Burke et al. 
[2^ predicted resonance II to occur at 356 ± 3 G which 
is 26 G lower than the observed value [l^. Taking into 
account that the experimental uncertainty in the posi- 
tion of the broad resonance II may be up to a few Gauss 



(it is not given in Ref. [la]), the agreement between our 
calculations and experiment can be considered very good. 

The resonances shown in Fig. [T] occur when the atoms 
are in weak field-seeking, magnetically trappable states. 
Optical traps, relying on the AC Stark shift, can provide 
confinement for arbitrary spin mixtures, and Feshbach 
resonances in high magnetic field-seeking states of *^Rb 
were recently observed at magnetic fields as small as a 
few Gauss [sllll. 

Figure [2] shows the s-wave Feshbach resonances in col- 
lisions of *^Rb in the 1) state with ^^Rb in three 
different hyperfine states: |2,0), |2, -1), and |2,1). All 
the resonances occur at magnetic fields below 10 G, and 
have widths of order 1-2 G (see Table |T|. They are simi- 
lar to the 9 G resonance found in a gas of pure *^Rb [s^l . 
We verfied that no other Feshbach resonances of similar 
widths occur at magnetic fields below 110 G. The reso- 
nances shown in Fig. [5] can be used to generate two-body 



spin entanglement in optical lattices [3l[, and to create 
heteronuclear molecules from atoms in different hyper- 
fine states. Finally, since these resonances occur at very 
low magnetic fields, they are technically much easier to 
access and they are possible candidates for the realiza- 
tion of RF induced Feshbach resonances using large RF 
or microwave fields generated near the surface of an atom 
chip |3a]. 



III. FESHBACH RESONANCES IN "LI- 
SYSTEM 
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The Li-Rb system is extremely important both from 
the standpoint of ultracold atomic gases and of ultracold 
molecular gases. Not only are ^Li and ^^Rb in widespread 
use for studies of ultracold fermionic and bosonic atomic 
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gases, but LiRb has a relatively large dipole moment. It 
is therefore an important candidate for the study of ul- 
tracold polar molecules and for the experimental study 
of electric- field-induced Feshbach resonances [265 . There- 
fore, understanding the low temperature coUisional prop- 
erties of this mixture is of great importance. 

In 2005, Silber et al. created a quantum degenerate 
Bose-Fermi mixture of ^Li and ^^Rb atoms in a rnag- 
netic trap with rubidium serving as the refrigerant [33| . 
This experiment revealed the challenges of this approach 
to cooling lithium due to small magnitude of the inter- 
species scattering length at low magnetic fields. Sub- 
sequently, inter-species Feshbach resonances in this sys- 
tem were found . Feshbach resonances may provide a 
means to enhance cooling in this mixture. Inter-species 
resonances also provide a way to tune the interactions 
in this Bose-Fermi mixture and may allow for the study 
of boson- mediated BCS pairing [34|. In addition, Fes- 
hbach resonances may offer an efficient way of forming 
loosely bound LiRb dimers, which can then be trans- 
ferred from the excited vibrational state near threshold to 
the ground vibrational state pol | . In deeply bound vibra- 
tional states, the LiRb dimer has a large electric dipole 
moment (of up to 4.2 Debye) [11], and an ensemble of 
these molecules, polarized by an external electric field, 
will interact strongly via the dipole-dipole interaction 
which is both long range and anisotropic. Such dipolar 
systems are predicted to exhibit a wide variety of novel 
phenomena [35] including superfluid, supersolid, Mott in- 
sulator, checkerboard, striped, and collapse phases for 
dipolar Bosonic gases [H, [s^] as well as novel superfluid 
phases of dipolar Fermi gases [H, [3§| and Luttinger liq- 
uid behavior in one dimensional traps (40| . 

Measurements of cross-thermalization in magnetically 
trapped ^Li-^^Rb and ''Li-^^Rb mixtures indicate that 
the interspecies triplet scattering lengths are lafrfjietl = 

20tl as m and \aJL^\ = 59t\l ob [H. In addition, 
two heteronuclear Feshbach resonances were recently ob- 
served [l^. The signs of the triplet scattering lengths 
and the location of Feshbach resonances in other atomic 
states, however, remains to be determined. The sign of 
the scattering length is particularly important since it de- 
termines the global stability of this mixture. Combining 
the experimental results we produce a new set of accurate 
LiRb interaction potentials which fully characterize the 
Li-Rb scattering properties in any combination of spin 
states and indicate that the sign of the ^Li-^^Rb triplet 
scattering length must be negative. Using these poten- 
tials, we also predict the location and widths of all the 
Feshbach resonances below 2 kG for all Li-Rb spin com- 
binations where '^^Rb is in the lower hyperfine manifold. 

A. Interatomic Potentials 

Our starting point for this work is to model the triplet 
a'^S and singlet X^E interaction potentials of LiRb by 
an analytical function of the form originally proposed by 



Degli-Esposti and Werner (i^ 

V{R) = G(i?)e-"(«-«=) - T{R) ^ ^ (6) 

1=3 

with G(i?) = ELoff'^' and T(i?) = i(l-htanh(l-hi-R)). 
The potential parameters were determined by varying 
this function to reproduce the overall shape and approx- 
imate number of bound states for the LiRb dimer ex- 
pected from the ab initio calculations [l^. The long 
range behavior is adjusted to match the van der Waals co- 
efficient Ce = 2545 Eh a% (where = 4.35974 x IQ-i^ J) 
determined by Derevianko et al. 24]. 

The amplitude and sign of the pure triplet and sin- 
glet s-wave scattering lengths as well as the positions 
of the Feshbach resonances are almost completely deter- 
mined by the location of the least bound states of the 
potentials Q. Since the long range behavior of the po- 
tentials has been accurately determined, the potentials 
can only be refined by making small adjustments to the 
short range repulsive wall while keeping the long range 
behavior fixed. This fine tuning is done to simultaneously 
reproduce the experimentally measured triplet scattering 
length and Feshbach resonance locations and widths. 

Since the full coupled channel calculation is computa- 
tionally intensive, iteratively finding the proper modifi- 
cation of the model potentials to reproduce the experi- 
mentally observed resonances can be a lengthy process. 
To simplify this search and to gain insight into the scat- 
tering properties of the Li-Rb system, we have employed 
the asymptotic bound state model (ABM) [IE O ISl 
to first determine the energies of the last bound states 
of the triplet and singlet potentials consistent with the 
experimentally observed location of the Li-Rb Feshbach 
resonances. We then tune the potential curves to repro- 
duce these bound state energies. Final refinement of the 
potentials is done to reproduce the exact location and 
widths of the observed Feshbach resonances. 



B. Results 

The ABM model is utilized and described elsewhere 
[TtI . and here we provide a brief summary of the model 
and the details of our use of it to facilitate the search for 
the correct interaction potentials. Feshbach resonances 
occur when a bound molecular state crosses the energy 
of the colliding atoms (open channel) at the dissociation 
threshold. The ABM model is therefore used to compute 
the energies of the molecular states (closest to the thresh- 
old) as a function of the magnetic field and to locate the 
crossings which result in Feshbach resonances. Since the 
Hamiltonian (in the absence of magnetic dipole-dipole 
coupling) preserves the total spin angular momentum, 
only states with the same Mp = Ms + Mj^ + Mj^ values 
as the initial state are considered. The major simplify- 
ing assumption in the ABM model is that the coupling 




Figure 3: Molecular bound state energies versus magnetic 
field computed with the asymptotic bound state (ABM) 
model. The threshold for the ||,^)6Li (8> |1, l)87pj^b collision 
channel is shown by the solid line while the dashed (dot- 
ted) lines indicate the s-wave (p-wave) states. These molec- 
ular state energies were computed given the least bound 
states, Eg, of the optimal singlet and triplet potentials of 
So"" = -0.106 cm-\ Eir^ = -0.0870 cm'^ , and = 
—0.137 cm~^, E{^^ = —0.116 cm~^. The predicted resonance 
locations are close to the actual locations determined by the 
full coupled channel calculation and are indicated by the solid 
dots. Near 890 G (A) the threshold crosses a p-wave molecu- 
lar state and the corresponding p-wave elastic scattering cross 
section shown in Fig.[S]is observed to rapidly diverge and then 
return to the background level. Likewise near 1070 G (B) the 
threshold crosses both a p-wave and s-wave molecular state 
and both the s-wave and p-wave elastic scattering cross sec- 
tions are affected. Finally, near 1300 G (C) a second s-wave 
induced Feshbach resonance occurs. 



between the channels (provided by the hyperfine interac- 
tion Vhf ) is small enough that the two-body bound states 
1 4"') can be represented to first order by uncoupled or- 
bital and spin states of the form j'l'') — \tpg) (g) |a) where 
\4>g) is the last bound state of either the pure singlet 
{S = 0) or triplet {S = 1) potential Vs{R). We there- 
fore replace the spatial part of the Hamiltonian ([T]) with 
Hid = J2s I ^sEg where Eg are the energies of these 
last bound states. In addition, we assume that the over- 
lap of the singlet and triplet wave functions of the same 
orbital angular momentum I is (V'olV'i) = !> the cou- 
pled bound state energies are found by diagonalizing the 
simplified Hamiltonian (iTj . 

Since the longrange part of the potential is known, the 
energies Eg^ of the higher I > states are uniquely de- 
termined by the I — triplet i?tripiot = E^ and singlet 
-Esingiet = -Eq energies, i^tripiet and i?singiet are therefore 
the only two free parameters in this model and they are 
adjusted until the threshold channel crosses the energy of 
the molecular states at the positions corresponding to the 
locations of the experimentally measured Feshbach reso- 



FIG. 4: The locus of points in the (-Esingict, ^'triplet) parameter 
space where an s-wave resonance occurs at one of the two ex- 
perimentally determined locations 882.02 G (grey/green) or 
1066.92 G (dark/red) for atoms in the |i, |)6Li (g) |1, l)87i^i, 
state. The dotted lines indicate the approximate values for 
-Esingiet and -Etripiot beyond which a new bound state enters 
the potential at zero energy. There are four regions (I-IV) 
indicated on the plot where an s-wave resonance occurs si- 
multaneously at 882.02 and at 1066.92 G. Region V indicates 
a range of values for which an s-wave resonance occurs at 
1066.92 while a p-wave resonance (not represented in this 
plot) occurs at 882 G. For each of these five candidate re- 
gions, the character of the predicted elastic cross sections as 
a function of magnetic field was studied and the results of this 
analysis are discussed in the text. 



nances. Figure [3] shows both the s and p-wave molecular 
bound state energies versus magnetic field for all states 
with Mp — 3/2 computed within the asymptotic bound 
state model. Although the ABM model cannot predict 
the exact location of the Feshbach resonances, it does 
predict reliably the energies of the molecular channels 
in regions far from the crossings. Therefore in the limit 
that the effect of the inter-state couplings on the energy 
is negligibly small, it provides an excellent estimate of 
the position of the Feshbach resonances. 

In order to utilize this model, we first employ it to 
determine the locus of points in the (iS'smgiet, ^'triplet) pa- 
rameter space where an s-wave resonance occurs at one 
of the two experimentally determined locations 882.02 G 
or 1066.92 G for atoms in the |i, ^)6u (8 |1, 1)^7 state. 
These points are plotted in Fig. 01 and four regions (I-IV) 
indicated on the plot are found for which an s-wave reso- 
nance occurs simultaneously at 882.02 and at 1066.92 G. 
Region V indicates a range of values (i?singict, -Etripict) 
for which an s-wave resonance occurs at 1066.92 G and 
a p-wave resonance occurs at 882 G. For each of the five 
candidate regions, the corresponding potentials were gen- 
erated and the predicted elastic scattering cross sections 
as a function of magnetic field were computed using the 
full coupled channel calculation. In addition, the corre- 
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spending triplet scattering lengths were also computed. 
Each of the four purely s-wave cases were ruled out based 
on a variety of reasons. In region I, the lower resonance 
at 882 G is predicted to be a factor of 10 larger in width 
than the upper resonance at 1067 G in violation of the 
experimentally measured widths of 1.27 G and 10.62 G 
respectively [la]. In region II, the relative widths of the 
resonances are (as in region I) incorrect and at these val- 
ues for (i^singiet, -^triplet) there would have been three ad- 
ditional and wide (> 5 G) s-wave resonances below 200 G 
which were not observed in the experiment. While in re- 
gions III and IV the ordering of the resonances is consis- 
tent with the experimental measurements (the upper res- 
onance is larger than the lower resonance), in region III 
there is an additional wide (>10G) s-wave resonance be- 
low 200 G not observed in the experiment, and in region 
IV there is an additional s-wave resonance at approxi- 
mately 960 G (>1G) in between the two observed reso- 
nances. In addition, the triplet scattering length for re- 
gions III and IV corresponding to i?tripict — —0.377 cm^^ 



IS a. 



6,87 
triplet 



iplct 

105 Ob. This value is in disagreement with 



the experimentally determined value from measurements 
of the cross-thermalization in magnetically trapped ^Li- 
^^Rb mixtures which indicate that the interspecies triplet 
scattering length is |a^;fJ[j,J = 20l:g as [111. 

In order to verify the robustness of these findings, we 
constructed a pair of Lennard- Jones potentials {V{R) = 
Ci2/R^^ — Cq/R^) with the same Cq coefficient, roughly 
the same number of bound states, and the same least 
bound state energies (£^singiet, ^^tripiet) as the fitted po- 
tential in each of the regions considered. Using these 
potentials, we verified that the Feshbach resonance loca- 
tions and scattering lengths are essentially the same as 
for the fitted potentials and are insensitive to the short 
range details of the potentials. This check provides an 
important verification of our characterization of the four 
purely s-wave candidate regions. The conclusion is that 
the experimentally observed Feshbach resonances are in- 
consistent with pure s-wave resonances, and we must con- 
sider the possibility that at least one of the resonances 
originates from a p-wave molecular state. 

Region V in Fig. 3] represents the only location in the 
(i^singiet, -^triplet) parameter space for which only one s- 
wave resonance occurs below 1.2 kG (at 1067 G) and a 
p-wave resonance occurs at 882 G. All other branches dis- 
played in Fig. [5] involve at least one additional s-wave res- 
onance occurring in a location where none was observed 
experimentally. Along the locus of (i?singiet, -^triplet) val- 
ues for which these two resonances occur at the correct 
locations, an additional p-wave resonance was found to 
occur somewhere between 1081 and 1024 G while the 
width of the s-wave resonance at 1065 G was found to 
vary from 5 G to 35 G. At the precise (i?singict, ^'triplet) 
values for which the second p-wave resonance was coin- 
cident with the s-wave resonance at 1065 G, the s-wave 
resonance width was AB = 11.53 G, consistent with the 

= 10.62 G, for 
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FIG. 5: Magnetic field dependence of the s-wave (upper 
panel) and p-wave (lower panel) elastic scattering cross sec- 
tions for atoms in the ||, 5)6]^; ® |1, i}»''B.h state. These re- 
sults are from the coupled channel calculations for a collision 
energy of 144 nK and using the optimal singlet and triplet po- 
tentials. Only the mi = contribution to the p-wave elastic 
scattering cross section is shown. Two s-wave resonances oc- 
cur at 1065 G and 1278 G while two p-wave resonances occur 
at 882 G and 1066 G. 



ture. For these optimal singlet and triplet potentials, 
the bound state energies are -Esingict = —0.106 cm^^ and 
-^triplet = -0.137 cm"^ 

Figure \5\ shows the result of the full coupled chan- 
nel calculation performed using the refined potentials. 
The elastic scattering cross-sections for the |^,5)6Li (8) 
|ljl)**''Rb state shows divergences at magnetic field val- 
ues at 1065 G and 882 G, in excellent agreement with 
the experimentally determined Feshbach resonance po- 
sitions. In addition, the triplet scattering length from 

6 87 

the fine tuned triplet potential was found to be aj^pjj,^ — 
— 19.8 Ob, also in excellent agreement with the experi- 
mentally determined value. For the reduced mass corre- 
sponding to a ^Li-®^Rb complex, the optimal fine tuned 



7,87 
triplet 



= 448 Ob, in disagree- 



triplct potential predicts a. 

ment with the experimental measurement of iLiti-ipj^ti — 
Close inspection of this potential reveals 
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41| 



experimentally measured value, ABexj 
the full width at half maximum for the trap loss fea- 



59tll OB 

that there is a bound state very close to the dissociation 
threshold for the ^Li-^^Rb triplet state. In the case, a 
small uncertainty in the exact location of this very weakly 
bound state (arising from uncertainties in the exact shape 
of the potential) translates into a very large uncertainty 
in the predicted triplet scattering length for ''Li-^^Rb 
mixtures. We verified the robustness of these results by 
changing the short range part of the potentials so that 
the number of bound states was different from the opti- 
mal potentials by more than 20% while still producing 
the same energy of the least bound states. As a result, 
the Feshbach resonance locations and scattering lengths 
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did not change significantly. In addition, we generated 
a set of Lennard-Jones potentials which reproduced the 
same least bound state energies and Feshbach resonance 
structure as the optimal fitted potentials. These poten- 
tials have a very different short range shape than the 
fitted potentials and they resulted in triplet scattering 
lengths of a^^^^^^^ = -22.6 ob and a[;fjj^^ = -333 ob con- 
firming that the determination for the ^Li--®^Rb triplet 
scattering length is very reliable (independent of the de- 
tails of the short range part of the potential) while the 
^Li-^^Rb triplet scattering length cannot be reliably pre- 
dicted given the proximity of a zero-energy resonance for 
this combination [i^ . We note that if the triplet scatter- 
ing length for the ^Li-*^Rb mixture is, in fact, negative, 
the experimental determination of its absolute magnitude 
can be complicated by the Ramsauer-Townsend effect in 
which the scattering cross section varies strongly with 
and may actually vanish at an experimentally relevant 
collision energy [44| . 

Using the refined potentials, the s and p wave scat- 
tering cross sections as a function of magnetic field were 
calculated for all spin combinations where ^^Rb is in the 
lower hyperfine manifold, and the location and widths of 
all resonances below 2 kG is summarized in Table HIl In 
experiments with ^Li-®^Rb mixtures, no Feshbach res- 
onances were observed below 1.2 kG for the ^)6u ^ 

|l,0)87Rb, I5, DsLi® |1, -l)smb, and |i,|)6Li«)|l,l)smb 
states. The results presented in Table|ll]are in agreement 
with the the last two of these observations but not the 
first. It is possible that because the resonances present 
in the |^, ^)<iLi ^ \^iO)^''Rh combination are very similar 
in location and width to those of the |i, ^}6Li'8' |1, l}s7Rb 
state that they may have been observed and erroneously 
concluded to arise from an impure state preparation. 

IV. CONCLUSIONS 

We have presented calculations of experimentally rele- 
vant magnetic Feshbach resonances for both ^^Rb-*^Rb 
and ^Li~®''Rb mixtures. Our results build on recent ex- 
perimental measurements of the triplet scattering lengths 
and Feshbach resonances in the Li-Rb system. We gen- 
erate a set of refined LiRb interaction potentials which 
reproduce the location and widths of the measured reso- 
nances with high precision and use these refined poten- 
tials to predict additional experimentally relevant res- 
onances for the ^Li-^^Rb mixture. These potentials 
indicate that the ^Li-^^Rb triplet scattering length is 
^trfpict ~ —19.8 flB, consistent with cross-thermalization 
measurements. We have verified that the predictions of 
these fine tuned potentials are robust in the sense that 
they only depend on the well-known longrange Cg co- 
efficient and are independent of both the details of the 
short range shape and the exact number of bound states 
of the interaction potentials. Finally, we have identified 
a zero-energy resonance generated by a bound or virtual 
state very near to the dissociation threshold in the ^Li- 



TABLE IL Position and width of '^Li-*''Rb Feshbach res- 
onances for magnetic fields below 2 kG determined from 
the coupled channel calculations. The experimentally mea- 
sured Feshbach resonances (and absence of resonances below 
1.2 kG) are also included for comparison. The experimentally 
determined width ABoxp is the full width at half maximum 
of the trap loss feature and, although related, it is not equiv- 
alent to AB (only defined for s-wave resonances). Several 
resonances were found which exhibited a suppressed oscilla- 
tion due to comparable coupling to inelastic channels [2^ and 
could not be assigned a width in the usual way (in accordance 
with Eq. O ■ In these cases, the maximum and minimum elas- 
tic scattering lengths of the oscillation were identified and the 
distance between them is indicated in parenthesis. 
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®'Rb triplet state. In the case of ^''^Rb-^'Rb collisions, 
using previously validated Rb2 interaction potentials, we 
have conducted a search for low field Feshbach resonances 
which would be candidates for the realization of RF in- 
duced Feshbach resonances [10]. The significance of this 
work is that low field resonances, below 10 G, do exist for 
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this mixture. These particular resonances are therefore 
very accessible and useful for chip based magnetic trap 
experiments. In addition, they may be eventually very 
useful for the generation of spin state entanglement in 
optical lattices or atom chip magnetic traps [45|. 
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